Spin Hall effect in iron-based superconducting materials: 
An effect of Dirac point 
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We have theoretically explored the intrinsic spin Hall effect in the iron-based superconductor 
family with a variety of materials. The study is motivated by an observation that, in addition to 
an appreciable spin-orbit coupling in the Fe 3d states, a character of the band structure in which 
Dirac cones appear below Ef may play a crucial role in producing a large spin Hall effect. The 
theoretical result does indeed predict a substantial spin Hall conductivity in the heavily hole-doped 
regime such as KFe2As2. The magnitude of the spin Hall effect has turned out to be comparable 
with that for Pt despite a relatively small spin-orbit coupling, which we identify to come form a huge 
contribution from the gap opening induced by the spin-orbit coupling at the Dirac point, which can 
become close to the Fermi energy for the heavy doping. 

PACS numbers: 72.25.Ba,74.25.Jb, 74.70.Xa, 75.47.-m 



There is a mounting interest in the spin Hall effect 
(SHE), where a transverse spin current (as opposed to 
charge current) is induced in an electric field [ij], which 
provides a most promising way for manipulating spin de- 
grees of freedom for spintronic devices based on nonmag- 
netic materials. After its theoretical prediction by Hirsch 
in 1999 [3|, which was originally proposed by D'yakonov 
and Perel' in 1971 [3[, investigations were initially fo- 
cused mostly in semiconductors [4]. More recently, atten- 
tions have been extended toward metallic systems |5l4l7{ , 
following experimental reports of a huge spin Hall con- 
ductivity (SHC) ~ 240&-in-icm-i in Pt ^^ that is 
about 10^ times larger than the reported values in semi- 
conductors y. A key question then is how we can iden- 
tify the materials that can exhibit large SHE. 

Although SHE is driven by the spin-orbit interaction 
(SOI), understanding its mechanism in metallic systems, 
especially intrinsic vs extrinsic mechanisms, remains an 
open issue [l|. While the intrinsic mechanism depends 
on the details of electronic structure, the extrinsic mech- 
anism is governed by impurities. Here we focus on the 
former, which has not onlysuccessfuUy accounted for the 
observed huge SHC in Pt 0, i, H El , but also predicted 
a varietyof materials, such as several Ad and 5d transition 
metals lfl,lll[, as good candidates for SHE. Interestingly, 
many of thesepredictions have been realized recently in 
experiments |16| . 

In this Letter we explore the possibility of SHE in the 
recently discovered iron-based superconducting materi- 
als, with their variety such as iron pnictides and chalco- 
genides. While the materials have become a class on its 
own of high- Tc superconductors after cuprates [18|, the 
reason why we have a fresh look at the materials from 
the viewpoint of SHE is the following. The iron-based 
superconductors are typically multiband systems, where 



various d orbitals are involved in the conduction bands 
|19| . The bands are entangled (namely, cross with each 
other), which gives rise to a situation where the orbital 
character changes as one walks along the Fermi surface. 
This is realized by an occurrence of "Dirac cones" as- 
sociated with the band crossing as an interplay between 
multiorbital character of bands and crystal symmetries 
[20|. We have also, for the iron compounds, an appre- 
ciable spin-orbit coupling (SOC) associated with the Fe 
3(i states, and a quasi-two-dimensional nature of disper- 
sion. Such features are expected to produce large SHE. 
For example, Kane and Mele have predicted a large and 
quantized SHC in graphene when the Dirac cone becomes 
massive due to an SOI-induced gap at the Dirac point [21 1 . 
However, due to extremely small gap due to a weak SOC 
(~ 4TOey)|22l for carbon 2p orbitals, its expermental re- 
alization remains challenging. By contrast, a relatively 
strong atomic SOC of order 72meV [23J for iron 3d or- 
bitals may make the SOC-induced gap at the Dirac points 
in the iron-based materials large enough for a realization 
of significant SHE. 



This is our reasoning, and, if the Dirac points in these 
materials lead to a large SHE, then, apart from being a 
good candidate for spintronics, the materials may also 
provide an avenue for exploring the Dirac physics in the 
context of SHE in metallic systems that has recently 
drawn a surge of interests in case of topological insu- 
lators [24]] . The present study does indeed predict a sub- 
stantially large SHC for heavily hole-doped system such 
as KFe2As2, which is contributed, as expected, by the 
electronic states in the vicinity of the SOC-induced gap 
at Dirac points that lie almost at the Fermi level for 
the heavily doped systems. We stress the Dirac-point- 
originated large SHE in the iron-based materials is dis- 
tinct from other systems such as Pt and 4d and 5d tran- 



sition metals, where substantially stronger SOC's govern 
the behavior of SHE. 

We consider the realistic band structure of various 
types of the iron compounds with two Fe atoms per unit 
cell by incorporating the SOC within an effective tight- 
binding (TB) Hamiltonian. The Wannier basis of the 
TB model has been constructed as follows. We first per- 
formed a density-functional calculation, where we used 
the exchange correlation functional proposed by Perdew 
et al. [25|, and the augmented plane wave and local 
orbital (APW-l-lo) method including the SOC as im- 
plemented in the WIEN2k codejll, ^. We then con- 
structed the TB model using the WIEN2 Wannier [H] and 
the wannier90 29]. In order to preserve the local sym- 
metry of the Wannier functions, we perform the so-called 
one-shot Wannierization [30|. Our approach provides a 
good description of several common electronic features 
in these materials, such as the pseudogap at the Fermi 
energy, an effective d-electron bandwidth of 4.5 — 5.0 eV, 
and the dominant orbital character of electronic states in- 
cluding the local d^z/dyz degeneracy, etc. These features 
themselves are almost independent of the SOC, as can 
be seen in Fig. [U and expected from the weak SOC for 
the Fe 3d electrons. However, we shall discuss later that 
some SOC-induced band features in fact play a crucial 
role in the SHE. 

SHE is investigated with the linear-response theorv[31j 
in the presence of local impurities that give rise to a finite 
residual resistivity. At T = 0, SHC consists of two parts, 
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represents the "Fermi-surface term" , while 
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represents the "Fermi-sea term" . Here the charge-current 
operator is given by Jp = ~edH /dk^ with — e(< 0) be- 
ing the electronic charge, and (^ — x,y, while the cr^-spin 
current operator is given by J^ — (— l/e){ JF, (T2}/2. 

The retarded (advanced) Green's function G^{G^) are 
given by G^/^(k,a;) = l/(w + ^i - H ± it) . Here H is 
the effective TB Hamiltonian, a 20 x 20 matrix spanned 
by the two Fe atoms per unit cell in the presence of SOC, 
and /i the chemical potential. F, the damping due to local 
impurities, is treated with the T-matrix approximation, 
F = (niinp/2«)[T(— iO) — T(-f iO)], where rtimp is the impu- 
rity concentration, and T(±iO) = I /[I — Ig{±iO)], with 
g{u}) = ^jj. G'(k, Lo) representing the T-matrix for a sin- 
gle impurity with an impurity potential /. For simplicity 
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FIG. 1: Band structures with (SOC) and without (nSOC) 
spin-orbit coupling are shown for KFe2As2. The symme- 
try points are: F = (0,0,0), X= (7r,0,0), Z= (7r,7r,0), 
T= (0, 0,7r), P= (7r,0,7r), and T'= (7r,7r,7r). In our con- 
vention, origin of energy is taken to be the chemical potential 
/I, and x~ and y~ are directed along the nearest neighbor Fe 
atoms. 



we consider a constant and orbital-diagonal impurity po- 
tential of strength /, i.e., Ia,f} = I5a,p- 

In the quantitative investigations that follow, our re- 
sults for SHC are expressed in unit of |e|/27ra, where 
a is the interlayer spacing. For a = 6A, \e\/2-Ka ^ 
645/ie~^r2~^cm~^. Unless otherwise mentioned, we con- 
sider / = 6 eV, and riimp = 0.01. 

Figure [2] shows the results for SHC for various materi- 
als, KFe2As2, FeSe, and LaFeAsOo.gFo.i, obtained with 
the ab initio band structures. The dominant contribu- 
tion is found to arise from the Fermi-surface part, i.e., 
a^y ~ a^]j, for all the materials considered. We can im- 
mediately notice that a large SHC arises in KFe2As2, 
a heavily hole-doped system, while it becomes vanish- 
ingly small in undoped FeSe and weakly electron-doped 
LaFeAsOo.gFo.i. Interestingly, the large magnitudes of 
SHC ^ 1300he-^n-^cm-^ in KFe2As2 with a resid- 
ual resistivity of '^ lO/xficm is comparable to that (~ 
1000/ie^^ri~^cm~^) predicted for Pt with a substantially 
larger SOC in the same metallic regime, where the ob- 
served SHC ^ 240fte~^ri~^cm~^ is reproduced at higher 
resistivity '^ lOO/ificm [TJ, |lOJ, [ill . We shall trace the ori- 
gin of the large SHC in KFe2As2 in a special band feature 
later. 

Since the SHE is found to be significantly affected by 
the doping level, we next investigate the effect of car- 
rier doping by systematically varying the band filling 
over a range {n — 5.5 — 6.5) that goes from electron 
doped (n > 6) to hole doped (n < 6) sides as relevant 
to the family of these materials. To single out the ef- 
fect of doping we adopt a rigid-band calculation using 
the band structure of BaFe2As2 that is known to cover 
both electron doping as in Ba(Fei_a;Coa;)2As2 and hole 
doping as in Bai_a;K3:Fe2As2. Here we only consider the 
Fermi-surface part, which is found to be dominant as dis- 
cussed above. The result, a curve in Fig. [51 exhibits a 
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FIG. 2: Results for the spin Hall conductivity for various 
systems (symbols). The curve is the doping dependence in the 
rigid-band approximation for the band structure of BaFe2As2. 
Inset shows the dependence of the result on the strength of 
impurity potential. 



large SHC in the heavily hole-doped regime (n ~ 5.5), 
while the magnitude is small in the weakly hole-doped 
and electron-doped regimes. The overall filling depen- 
dence agrees with the results for KFe2As2 {n — 5.5), 
FeSe (n — 6.0), and LaFeOo.gFo.i {n — 6.1), and high- 
lights the presence of a common feature. Small deviations 
from the rigid-band result should be due to deformations 
in the electronic structures, such as the stronger three 
dimensionality in BaFe2As2. 

As for the dependence on the impurity potential 
strength, we have checked, as shown in the inset of Fig. 
[21 that, when we vary the strength over (/ ~ 2 — 20) eV, 
the SHC, including its filling dependence, is rather in- 
sensitive to the strength of impurity potential [32|. This 
hallmarks an intrinsic nature of SHE in these materials. 

Let us now focus on the case of heavily hole- 
doped KFe2As2 that exhibits the large SHC. In or- 
der to identify its origin we explore if certain special 
band features around the Fermi energy are responsi- 
ble. For this purpose we look into the momentum- 
dependent contribution of electronic states on the 
Fermi surface to the dominating Fermi-surface part in 
SHC, cr^yOi.) = 8 X]fc'=±fc,,±fc„,±fe, Tr[(J^G' Jy G - 
Jy G'^ J^ G^) /2]k\ij=o where we are averaging over xy 
and yx components etc, with 2ilv Sk'''^j/(k) providing 



the net contribution of the Fermi surface part a^]^ [7| . 

We have only to consider a part of the Brillouin zone, 
as shown by the the dotted area in Fig. [3] (a), while 
the other parts can be deduced from the symmetry. For 
clarity we consider the kz = n plane, as other planes 
give qualitatively similar contributions due to quasi-2D 
nature of the Fermi surface that is also apparent from 
the almost flat dispersion near the Fermi level along F — 
T direction (Fig. [T]). Figure |3] (b) shows that a huge 
contribution to SHC comes from a small region on the 
hole pocket near P. On the other hand, contributions 
from the rest of the Fermi surface are not only small but 
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FIG. 3: (a) Fermi surface for KFe2As2 shown on the k^ = 
TT plane, (b) k-dependent contribution from the states on 
the Fermi surface (red lines) to SHC (green), displayed for 
a part of Fermi surface (dotted region in (a)). Contribution 
from the other parts of the Fermi surface can be deduced 
from the symmetry, (c) Gapped Dirac cone dispersion near 
P. Electronic states make the huge contribution to SHC (b) 
in the close vicinity of the gap. 



also tend to cancel with each other. Therefore, origin of 
the large SHC can be attributed to the huge contribution 
of electronic states on the hole pocket near P. 

We can actually identify, in Figures |3] (c) and SI that 
these electronic states lie in the close vicinity of a small 
SOC-induced gap in the Dirac cone, which is tilted. This 
is also the case with our analysis for the filling depen- 
dence in BaFe2As2 with a rigid band, where we find a 
similar origin behind the large SHC in the heavily hole- 
doped regime {n ~ 5.5), in which the Fermi level ap- 
proaches the Dirac point. Concomitantly, such Dirac- 
like points near the Fermi level are absent in FeSe and 
LaFeAsOo.gFo.i. Such Dirac-like points arise due to 
crossing of two bands that respectively have almost lin- 
ear dispersions and different orbital characters as a con- 
sequence of crystal symmetries in the iron-based mate- 
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FIG. 4: The band dispersions for the electronic states in 
KFe2As2 without SOC (a) and with SOC (b). The Dirac 
point, around which dominant contributions to SHC arise as 
shown in Fig. [3l is indicated with arrows. 



rials, as also discussed recently in detail by Anderson 
and Boeri[2(|. In the present case, as shown in Fig. |4j 
while one of the crossing bands has dominantly yz or- 
bital character, the other band has xy character. The 
band crossing is a consequence of the opposite parities of 
the two orbitals with respect to reflection to xz mirror 
plane. In the presence of an SOI, the gap opens at the 
Dirac point because the two orbitals are then coupled by 
the y— component of the orbital angular momentum ly] 
{yz\iy\xy) = -i. 

Thus the large SHC when the Fermi level lies close 
to the SOC-induced gap at the Dirac point renders 
the iron-based systems a Dirac electron systems, as in 
graphene. 2l| The difference here is, unlike in graphene, 
the gap can be large enough to make the SHE experi- 
mentally accessible due to relatively strong strength of 
SOC for 3d electrons. We briefly comment on the role of 
current vertex corrections due to local impurity poten- 
tials. In the present investigation, as we consider only Fe 
3d electrons, they vanish identically as a consequence of 
the even parity of d-electrons with respect to k — ;• — k 
[lO|. However, in the presence of p electrons with an odd 
parity that come from the other elements such as pnicto- 
gens and chalcogens, the current vertex corrections may 
be finite, as indeed found in graphene[33J. Nevertheless, 
we expect that the present result would not be modified 
significantly as the dominant contribution to electronic 
states around the Fermi level comes from the Fe 3d or- 



bitals. 

In summary, we have theoretically explored the pos- 
sibility of spin Hall effect in a variety of the recently 
discovered iron-based superconducting materials. We re- 
veal that a substantially large spin Hall current arises 
in heavily hole-doped systems such as KFe2As2, whose 
magnitude is comparable with that in Ft. The large spin 
Hall current is found to come from a huge contributions 
from the electronic states in the vicinity of the spin-orbit 
coupling induced gap at the Dirac points that is made to 
lie close to the Fermi level in the heavily doped case. 
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